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Abstract 

> , 

o : 

' We describe a new wave mode similar to the acoustic wave in which both density and velocity fluctuate. 

Unlike the acoustic wave in which the underlying distribution is Maxwellian, this new wave mode occurs when 
the underlying distribution is a suprathermal kappa function and involves fluctuations in the power law index, 
K. This wave mode always propagates faster than the acoustic wave with an equivalent effective temperature 
I ' and becomes the acoustic wave in the Maxwellian limit as k — > oo. 

o ■ 

, 1. Introduction 

' Distributions in space plasmas are generally described by a drifting Maxwellian, 

= ^^^expl- *-^ ""^ }, (1) 

in one-dimension. This type of distribution supports wave modes in which the various physical parameters: n, 
1 the density, uq, the thermal speed, and u, the convection speed, can all oscillate. Perhaps the simplest of these 
^ " is the acoustic or sound wave. 

■ However, space plasma distribution functions are not Maxwellian [e.g. references in Collier, 1993], but 

H : 

' suprathermal, generally well-described by a drifting kappa function which in one dimension is [Summers and 
Thorne, 1991]: 

WoVTT't'^^ r(K- i) [1 + (v - U) / KU^] 

Here, the parameters n and u are identical to the above, loq is a parameter analogous to that in equation (1), 
and K is an additional parameter describing the prominence of the suprathermal tail, the lower the value of k, 
the more pronounced the tail. In the limit as k ^ oo, equation (2) becomes equation (1). 

Although many view (2) as a non-equilibrium distribution, nature contradicts this notion. Suprathermal 
distributions well-described by kappa functions occur almost universally in space plasmas [e.g. Summers and 
Thorne, 1991; Collier, 1999] and theoretical considerations [Summers and Thorne, 1991; Leubner and Voros, 
2005; Leubner, 2004a; Treumann et al., 2004], for example based on the Tsallis entropy [Tsallis, 1995], have 
aided in understanding their preeminence. In this manuscript, we shall adopt the point-of-view that the kappa 
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function, rather than the Maxwellian, is the equilibrium distribution. Of course, we recover the Maxwellian 
results when k — > oo. 

We have come to expect wave modes to involve the parameters of physical relevance identifiable in the 
Maxwellian distribution, namely the moments density, temperature (or equivalently uq), and flow (or u). How- 
ever, there exists no physical reason why wave-like behavior should be limited solely to moments of the dis- 
tribution. Because the kappa function, in addition to containing the associated aforementioned Maxwellian 
parameters, has the extra parameter k, it seems reasonable, then, indeed expected, that for suprathermal distri- 
butions there should also be wave modes in which k participates. Here we provide an example of such a wave 
mode in a one- dimensional neutral fluid based on the ion acoustic wave. 

Ion acoustic waves play an important role in the interaction between bodies moving supersonically through 
a plasma and the ambient plasma itself. In this case, they create a wake such as in the region behind Venus, 
the moon, lo and Titan. In the wake, converging beams of flowing plasma collide, exciting ion acoustic waves 
[e.g. Samir et al., 1983]. Broadband electrostatic noise that could plausibly be identified with ion acoustic wave 
activity has been observed in the polar cusp [Gurnett and Frank, 1978], the distant magnetotail [Gurnett et al., 
1976] and along auroral field lines [Gurnett and Frank, 1977]. 



2. The Collision Operator 

Although the presence of suprathermal tails may be interpreted as evidence of non-equilibrium distribu- 
tions, the near-universal occurrence of these distributions in space plasmas [e.g. Summers and Thorne, 1991; 
Collier, 1999] argues in favor of a formalism which treats these distributions as a physical equilibrium (but not 
thermodynamic equilibrium) state rather than as an anomaly. 

To this end, we consider the one-dimensional force-free Boltzmann equation given by 



dt dx dt 



(3) 



where / is the distribution function, v is the particle velocity, x is position, t is time, and df/dt\^ is the collision 
operator. 

The standard approach is to assure that collisions reduce the distribution to a Maxwellian consistent with 
the local thermodynamic properties. However, in the case considered here, because the equilibrium distribution 
becomes the K-function, the effect of the collision operator (under the appropriate constraints to be discussed 
shortly) must be consistent with this distribution. 

Some properties of the collision operator remain unchanged in spite of the K-function assumption. Because 
collision operators are "local" in configuration space, if two particles collide, they make a jump in velocity space 
but remain at approximately the same physical location before and after the collision. Thus, 



a/ 

dt 



dv = 0, (4) 
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because at a given position even though the particles there shift around in velocity space due to collisions, the 
local number density must not change. Additionally, because momentum must be conserved in any collision 



f 

J — c 



5/ 
dt 



vdv = 0. (5) 



The standard treatment also assumes conservation of energy. However, there exists no reason to expect 
energy to be conserved in whatever process establishes the equilibrium kappa distribution. It can be shown 
that the assumption of a constant order of magnitude of energy < ln{H} >, where H = ^ + ^mv"^ with $ the 
(positive) potential, m the mass, and v the velocity, leads to an equilibrium distribution that is a kappa function 
[Collier, 2004]. Thus, the collision operator obeys 



r 

J — ( 



dv = 0. (6) 



One final note: this approach assumes that collisions arc dominant and that the collision frequency is 
the highest frequency for the phenomenon under consideration. Of course, this is an invalid assumption for 
collisionless space plasmas. Nevertheless, we justify this assumption on the basis that this is the standard 
approach for deriving the fluid equations, which work far better than they should for space physics applications. 



3. The Equation Set 

The first two equations with which we will work are 

dn d , , ^ 

du dP d . K 1 2\ tr.\ 

Equations (7) and (8) are simply the zeroth and flrst moments of equation (3) which represent physically mass 
conservation and momentum. The time derivative on the left of equation (8) is the convective derivative. In 

equation (8), wc have substituted for the pressure, P, the product of the density and temperature. Here we 
introduce an artificial "core" temperature, Tcore^ defined by 



-gTcore = 3 ' -WWo- (9) 



K- 2 ■"2 

The core temperature is simply the corresponding Maxwellian temperature for k — > oo and is not meant to imply 
that the distribution has two temperatures since the temperature, a moment of the distribution function, can be 
calculated uniquely for any distribution [Leubner, 2004b]. Note that the factor k — | in the denominator appears 
in both the one-dimensional and three-dimensional cases. 

The justification for defining the variable Tcore is that the expression for the entropy of the K-function: 

Tcire 

— — — ^id(k) = constant, (10) 
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where 



^'^^''^ = ''r[l + 1)'^ exp|«[V'(«) - V(« - 1/2)] - 1/2|, 



(11) 



and tp is the psi or digamma function equal to the derivative of the logarithm of the gamma function [Abramowitz 
and Stegun, 1972], will be used to close the fluid equations. Because this expression involves separate factors 
depending only on uq and k respectively, defining the additional variable Tcore will prove convenient for the 
following discussion. The function is shown in Figure 1 along with the three dimensional case for comparison 

[Collier, 1995]. 

Equation (10) results from the ln{/}/n moment of equation (3), noting that the equilibrium distribution 
is a kappa function, that is using equation (6), and is physically conservation of entropy [Collier, 1995]. To 
achieve equation (10), we consider f{x,v,t) = n{x,t) ■ g{x,v,t) where g{x,v,t) is the normalized velocity space 
distribution, that is, 

/oo 
g{x,v,t) dv = l. (12) 
-OO 

Taking the hi{ng}/n moment of the right hand side of equation (3), the colhsion operator, gives us 



dl 
dt 



1 , r . , I f°° df 

-\n{ng}dv = - — 
n n dt 



H9}dv, (13) 



because of conservation of particles as expressed by equation (4). 
The first term on the left hand side may be expanded as 



f°° df I d f f dn f dn 

/ ■ -'^'^{'"'9}dv = —{-\n{ng})dv+ {—) — \n{n g} dv - {—)^ 



dv, (14) 



where we have used the fact that the integral of the time derivative of g is zero because g is normalized to unity, 
equation (12). 

Likewise, the second term on the left hand side may be expanded as 



/■°° 5/ 1 , r . , r d ,fv ^ , , Jv.dn^ , , , 

J-oo dx 



,v f .dn , f°° do 
y dx 

In the sum of equations (14) and (15), the last term in equation (14) and the two last terms in equation (15) 
combine to yield 



dn dn du 
dt dx dx 



0, (16) 



where the terms in the square brackets sum to zero by mass conservations, equation (7). 
Consequently, we can write the left hand side of the ln{/}/n moment of equation (3) 

P 5/ 1 1 r ^ , r df 1^ , , ^ 
/ • -M"ffl«'*'+ / V- ln{ng}dv = 

J-oo dt n dx n 

/CO o r pOO f f) fi pOO r pOO r o 
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where the third term on the right hand side arises from the assumption that / is symmetric around u so that 

/ (v-u) -ln{ng}dv = 0, (18) 

which is equivalent to asserting that there is no heat flux. Equation (17), then, may be re-written as 

f°° df 1 f°° df 1 

j^[J_jHn9}dv] +J_-.Hn9}[^+n-+u-]dv, (19) 

where d/dt is the convective derivative, and we have again used equation (18). 

Noting that the terms in the square bracket in the second integral on the right hand side of equation (19) 
sum to zero by mass conservation, equation (16), we arrive at 



^r^{9}dv, (20) 



where we have used the right hand side of equation (13). 

As discussed in Section (2), we assume that the collision operator effects an equilibrium kappa distribution 
in accordance with Collier [2004] so that 

g{v)<x{^+^mv^) . (21) 
Under this condition, the right hand side of equation (20) is zero, so that the integral 

/oo 
g ln{n g} dv = constant. (22) 
-co 

Using the g defined by the one-dimensional kappa function, equation (2), noting that 



d_ 

da 



using the identity 



and the definition of 



^=T{zMz), (25) 



one retrieves equation (10) which is the entropy of the one-dimensional kappa function. Note that using the 
large z asymptotic expansions 

r(^) « \/2^exp{-^}2^-V^ (26) 

V(^) « Hz} - 1, (27) 
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one can show that 



6 



lim aD(«) = l, (28) 
so that equation (10) becomes the standard expression for conservation of entropy in a Maxwellian gas. 

4. The Dispersion Relation: Standard- Acoustic Mode 

Summarizing Section 3, the set of equations with which we will work are 

dn d r T ^ ,. N 

^ + ^[H=0, (29) 

rpl/2 

——^idM = const, (31) 
n 

and we will apply the standard approach of linearizing these equations assuming sinusoidal variation of the 
relevant parameters. 

Equation (29), mass conservation, requires that both density n and velocity u oscillate. So, we consider 
the standard equilibrium in which density oscillates around some equilibrium Uq and ignore any equilibrium 
convection. However, both T^ore o,nd k need not oscillate simultaneously so that this set of equations supports two 
distinct modes. We will consider first the mode in which core temperature Tcore oscillates but the suprathermal 
tail K does not so that 

n = no + 6n exTp{i{kx — Lot)} , (32) 
u = 6u exp{i{kx — cut)} , (33) 
+ STcore exp{i(fca; -ujt)], (34) 

where the ^'s correspond to first order perturbations. 

Expanding equations (29) through (31) to first order using equations (32) through (34) gives us 



io5n = uokSu, (35) 

k, (36) 



m K - I 



coreo 



1 STcore Sn 



2 Tqqj-qq 77-0 



(37) 



Combining equations (36) and (37) to eliminate 5Tcore/Tcoreo and then using equation (35) to eliminate lo5u 
eventually gives us 



(38) 
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which is simply the dispersion relation for the one-dimensional acoustic wave [e.g. Chen, 1984] except that, in 
this case, the underlying distribution is a kappa function with suprathermal tails that simply "go along for the 
ride," so to speak. Though, note that T in equation (38) above is the total temperature, not core temperature, 
so that the suprathermal character is still important in that sense. 



5. The Dispersion Relation: Kappa- Acoustic Mode 

Now we will make the assumption that only density, velocity and kappa oscillate according to 

n = uq + Sn exp{i(kx — cot)}, 

u = (5uexp{j(fca; — ojt)}, 

K = kq + (5Kexp{i(fca; — wt)}, 

where the 5s correspond to first order perturbations. 

Using equations (29), (30) and (31), we get to first order 



(39) 
(40) 
(41) 



uj5n — nokSu. 

e '-^core ^0 
UJOU = o- 

m Ko - I 

5n 



\5n 



Sn 



no ko{ko - i) 



H = — 

OK no 



Combining equations (43) and (44), we get 



jSu 



T 



m Ko — 5 

Combining equations (42) and (44), we get 

ivSu = 



dn 



ko{k,o - 1) 



5k. 



6k,. 



(42) 
(43) 

(44) 



(45) 



(46) 



k Ok 

Finally, combining equations (45) and (46), we arrive at the dispersion relation for the one-dimensional kappa- 
acoustic wave: 



/ST I 1 1 
^ V3~2'^„(^P_3) Ointe^W T 



(47) 



Note that w in equation (47) is never imaginary. 



6. Discussion 

By examining equation (47) and noting that the phase speed of a standard one-dimensional acoustic wave 
Cs with a temperature T is 
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we can define the ratio of the phase speed of the kappa-acoustic wave to the phase speed of a standard acoustic 
wave of the same temperature as 



^-\l3-2\,^,^_s^eM^- (49) 

By using the definition of £,id{k) given by equation (11) and noting that ^P{k + 1) — fp{K) = 1/k, we can show 
that 

Here tp' , the derivative of the gamma function, is referred to as the trigamma function [Abramowitz and Stegun, 
1972]. Figure (2) shows equation (49) plotted as a function of kappa. The phase velocity of the acoustic kappa 
wave is always greater than a standard acoustic wave of the same temperature and the two phase velocities 
are identical in the limit of large k, as can be shown by using the asymptotic expression for large k, t/j'in) 
1/k + 1/2k;^. The phase velocity of the kappa wave diverges to infinity as the value of kappa goes to 3/2. 

In addition to the obvious characteristic that in the kappa acoustic mode the suprathermal tail oscillates, 
note that for the kappa acoustic wave, the relative density fluctuation obeys 

6n dlnl^io} 



no Ok 

while for the standard-acoustic mode the relationship is 

Sn 1 STc, 



-5k. (51) 



(52) 



^0 2 ^coreo 

SO that in the kappa acoustic mode, the compressions in density are associated with decreases in k or more 
pronounced tails. This, then, is the equivalent of "heating" the distribution except that the energy goes into the 
tail rather than the core. However, the restoring force is different because the degree to which the tail can be 
tweaked is different than that of the core and is given by the equation of state, equation (10). 



7. Conclusions 

We have described a new wave mode which is analogous to the acoustic wave, but involves oscillations of the 
tail of suprathermal distributions. These waves always travel faster than sound waves with the same temperature 
and in the limit as k; — > oo become standard acoustic waves. 

The simple analysis presented here is based upon a one- dimensional collisionless neutral gas. In actuality, 
most physical contexts of interest are more complex, involving plasmas with electric and magnetic fields. However, 
the purpose of this paper is primarily to illustrate that investigators should consider the possibility and look for 
evidence in the data of wave modes involving the suprathermal tail, that is k. 

Recently, however, there has been some debate about low frequency (few Hz - kHz) broadband electrostatic 
emission representing ion acoustic wave activity versus solitary waves triggering a wide response in the waveform 
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capture. One way to resolve this issue, at least at the lower frequencies below about 10 Hz, is to take in addition 
to electric and magnetic field measurements, complete ion distribution functions at frequencies close to 10 Hz 
or a cadence of 100 ms. Ion acoustic waves would be identified by correlations between density and velocity, In 
addition, possible non-Maxwellian generalizations of the ion acoustic wave, which are the subject of this paper, 
could be observed with such instrumentation. 

Unfortunately, current magnetospheric particle instrumentation has never had a sufficient cadence to observe 
the new wave mode predicted in this paper. One possibility is that the high-time resolution plasma investigations 
on the Magnetospheric Multiscale (MMS) mission currently scheduled to launch in 2014 may be able to detect 
these types of waves. The MMS/SMART Fast Plasma Instrument (FPI) Dual Ion Spectrometers (DIS) will 
acquire complete distribution functions up to 40 keV every 37.5 ms (four energy sweeps every 150 ms). The 
orbit during the first of the three MMS phases finds the four spacecraft formation in a 1.2x12 Re orbit at 28° 
inclination which will permit high-time resolution studies of both magnetospheric and solar wind plasmas. 

Acknowledgments. Thanks to Bill Farrell, John Podesta, and John Sigwarth for helpful conversations. 
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Figure Captions 

Figure 1. The one and three-dimensional ^ functions which modify the standard Maxwellian expressions for 
conservation of entropy. 

Figure 2. The ratio of the phase velocity of the one-dimensional acoustic kappa wave normalized by the phase 
speed of a standard acoustic wave with the same temperature. 



